load belastingsbussen  %gegevens over de verschillende belastingsbussen laden
load problemen
kosten_mitigatiemiddelen;
sectoren_namen; %namen sectoren laden
[NUMERIC1,TXT,IndelingBussen] = xlsread('Verdeling_belastingsbussen.xls');
[r,k]=size(IndelingBussen);
kost_deze_bus_miti=0;
kost_deze_bus_PQ=0;
totale_kost_miti=0;
totale_kost_PQ=0;
for i=1:1:aantal_belastingsbussen
   nummer_belastingsbus = belastingsbussen(i,1);
   rijen=find(cell2mat(IndelingBussen(2:size(IndelingBussen,1),1))==nummer_belastingsbus)
   nummer_belastingsbusstr = int2str(nummer_belastingsbus);
   aantal_sectoren_op_deze_bus=size(rijen,1)
   for l=1:1:aantal_sectoren_op_deze_bus
       rij=rijen(l,1);
       naam_sector=char(IndelingBussen(rij+1,3));%Plus 1 want in de eerste rij van het excel sheet staan titels
       test_string=strcat('kosten_sector_',naam_sector);
       aandeel_deze_sector=cell2mat(IndelingBussen(rij+1,4))/100;
       temp=char(test_string);
       test=eval(temp);
       kost_deze_sector_miti=test(i,3)*aandeel_deze_sector;
       kost_deze_sector_PQ=test(i,2)*aandeel_deze_sector;
       kost_deze_bus_miti=kost_deze_bus_miti + kost_deze_sector_miti;
       kost_deze_bus_PQ=kost_deze_bus_PQ + kost_deze_sector_PQ;
   end
   totale_kost_miti=totale_kost_miti + kost_deze_bus_miti;
   totale_kost_PQ=totale_kost_PQ + kost_deze_bus_PQ;
   Kosten_Per_Belastingsbus(i,3) = kost_deze_bus_miti;
   Kosten_Per_Belastingsbus(i,2) = kost_deze_bus_PQ;
   Kosten_Per_Belastingsbus(i,1) = nummer_belastingsbus;
   baten = Kosten_Per_Belastingsbus(i,2) - Kosten_Per_Belastingsbus(i,3);
   ratio = baten/Kosten_Per_Belastingsbus(i,3);
   Kosten_Per_Belastingsbus(i,4)=ratio;
   kost_deze_bus_miti=0;
   kost_deze_bus_PQ=0;
end
%kost_deze_bus
Kosten_Per_Belastingsbus
save Kosten_Per_Belastingsbus;
